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We propose a physically motivated parametrization for the unpolarized generalized 
fT^ I parton distributions, H and E, valid at both zero and non-zero values of the skewness 

variable, ^. At — 0, H and E are determined using constraints from simultaneous 
fits of experimental data on both the nucleon elastic form factors and the deep in- 
elastic structure functions. Lattice calculations of the higher moments constrain the 
parametrization at ^ > 0. Our method provides a step towards a model independent 
l' extraction of generalized distributions from the data that is alternative to the mathe- 

O , matical ansatz of double distributions. 

1 Introduction 

>, 

, Generalized Parton Distributions (GPDs) parametrize the soft contributions in a variety of 

hard exclusive processes, from Deeply Virtual Compton Scattering (DVCS) to hard exclu- 
sive meson production (see [2j [3] for reviews). The conceptual idea behind their definition 
allows one to address a vast, previously inaccessible phenomenology, from the simultaneous 
description of hadronic structure in terms of transverse spatial and longitudinal momentum 
degrees of freedom [3] , to the the access to the description of angular momentum of partons 
, in nucleons and nuclei via Ji's sum rule [5]. 

At present, a central issue is the definition of a quantitative, reliable approach beyond 
the construction of GPDs from specific models and/or particular limiting cases, that can 
' incorporate new incoming experimental data in a variety of ranges of the scale Q^, and 

^ , the four- momentum transfer between the incoming and outgoing protons, A = (i, C)- The 

matching between measured quantities and Pcrturbative QCD (PQCD) based predictions for 
DVCS should proceed, owing to specific factorization theorems, similarly to the extraction 
of Parton Distribution Functions (PDFs) from deep inelastic scattering. A few important 
caveats are however present since GPDs describe amplitudes and are therefore more elusive 
observables in experimental measurements. Experiments delivering sufficiently accurate 
data have, in fact, just begun [6]. The comparison with experiment and the formulation 
of parametrizations necessarily encompasses, therefore, other strategies using additional 
constraints, other than from a direct comparison with the data. 

We propose a strategy using a combination of experimental data on nucleon form factors, 
PDFs, and lattice calculations of Mellin moments with n > 1. The latter, parametrized in 
terms of Generalized Form Factors (GFFs), were calculated by both the QCDSF j7j and 
LHPC [8] collaborations for both the unpolarized and polarized cases up to n = 3, therefore 
allowing to access the skewness dependence of GPDs. 



o 
p 

o 



DIS2007 



2 Generalized Parton Distributions from Lattice Moments 



GPDs can be extracted most cleanly from 
DVCS [5]. In this contribution we concen- 
trate on the unpolarized scattering GPDs, 
H, and from the vector (7^) and ten- 
sor {(T^iu) interactions, respectively. We 
adopt the following set of kinematical vari- 
ables: (C,^,i), where C = Q^/2{Pq) is the 
longitudinal momentum transfer between 
the initial and final protons (C ~ xsj in 
the asymptotic limit, with Bjorken xsj = 
Q'^/2Mv), X = {kq)/{Pq) is the momen- 
tum fraction relative to the initial proton 
carried by the struck parton, t ~ — A^, is 
the four-momentum transfer squared. X is 
not directly observable, it appears in the 
amplitude as an integration variable [3] . 
The need to deal with a more complicated 
phase space, in addition to the fact that 
DVCS interferes coherently with the Bethe- 
Heitler (BH) process, are in essence the rea- 
sons why it is more challenging to extract 
GPDs from experiment, wherefore guid- 
ance from phenomenologically motivated 
parametrizations becomes important. 
We first present a parametrization of H and 
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Figure 1: (color online) The dipole masses 
squared for n = 2, for the isovector magnetic 
(lower panel) and electric (upper panel) con- 
tributions obtained by performing fits to the 
lattice results of [7| . The value at the physical 
pion mass obtained from our fit is also shown 
(star). 



E in the flavor Non Singlet (NS) sector, 

valid in the X > C, region, obtained by extending our previous zero skewness form [9], 
through proper kinematical shifts: 



H{XX,t) = G\i^{XX,t)R{XX,t) 



(1) 



(a similar form was used for E{X,C,,t)), where R{XX,t) is a Regge motivated term de- 
scribing the low X and t behaviors, while G^.[^{XX,t), was obtained using a spectator 
model. 

In order to model the X < ( region, we observe that the higher moments of GPDs 
give C-dependent constraints, in addition to the ones from the nucleon form factors. The 
n = 1, 2, 3 moments of the NS combinations: i/"^'' ^ H" - H'^, and = E"" - E"^ are 

available from lattice QCD [71 [8] . They can be written in terms of the isovector components 
as: 



e: 



L—d 



= / dXX' 



11^''^ = J dXX"-^{E'' ~ E"^) = 



1 



1 



(2) 
(3) 



where the l.h.s. quantities are obtained from the lattice moments calculations, whereas 
{HXd(E))^ (-^JU(£;))" amenable to chiral extrapolations. We used lattice calculations 
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for the unpolarized GFFs obtained by the QCDSF collaboration using two flavors of 0{a)- 
improved dynamical fermions for several values of t in the interval < t < 5 GeV^, and 
covering a range of pion mass values, 171^^ > 500 MeV^. Similarly to previous evaluations 
[8] the GFFs for both H and E, display a dipole type behavior for all three n values, the 
value of the dipole mass increasing with n. We performed an extrapolation by extending to 
the n = 2, 3 moments a simple ansatz proposed in jlO] for the nucleon form factors that: 
i) uses the connection between the dipole mass and the nucleons radius; ii) introduces a 
modification of the non analytic terms in the standard chiral extrapolation that suppresses 
the contribution of chiral loops at large m^r. Despite its simplicity, the ansatz seems to 
reproduce both the lattice results trend at large m^r while satisfying the main physical 
criteria i) and ii). Our results for the dipole mass at n = 2 are shown in FigUJ 



3 Reconstruction from Bernstein Polynomials 



Similarly to the PDFs case [TT], with a 
finite number of moments in hand, one 
can use reconstruction methods attaining 
weighted averages of the GPDs, around av- 
erage ranges of X. The weights are provided 
by the complete set of Bernstein polynomi- 
als. 

In Figl2]we show H^^'^ reconstructed us- 
ing the available lattice moments. We per- 
formed the procedure in the X < C, region 
only using: 



C 



(4) 



Hk,niC,t)= J H{XX.t)huAXX)dX 


where: 6fe,„(X,C) = X^ (C-X)"-V (C- 



and we used subtracted mo- 
ments, defined as: 
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Figure 2: (color online) Comparison of i/" 
for different values of ^and —t = tmin = 
0.035,0.073,0.18,0.53 GeV^ (left panel), and 
— t = — 1 GcV^ (right panel), calculated using 
the procedure described in the text. 



where Hn are the Mellin moments, and 
H'{X,C,t) was obtained from Eq.((l|). For 
71 = 2, fc = 0,l,2, the reconstruction proce- 
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dure yields [T^ : 



H 22{C,X22) 




A 



32 



(6) 



where Xqi = 0.25, X02 = 0.5, X03 = 0.75, and ^10,^20,^30, ^32 are the GFFs from Ref.[7]. 

In conclusion, we provided a fully quantitative parametrization of the NS GPDs, valid in 
the region of Jefferson Lab experiments [3] that, differently from model calculations, and for 
the first time to our knowledge, makes use of experimental data in combination with lattice 
results. Given the paucity of current direct experimental measurements of GPDs, our goal is 
to provide more stringent, model independent predictions that will be useful both for model 
builders, in order to understand the dynamics of GPDs, and for the planning of future hard 
exclusive scattering experiments. 
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